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The nonlinear interaction between a pair of symmetric, oblique, and spatial instability 
modes is studied in the long-wave limit using asymptotic methods. The base flow is 
taken to be a supersonic mixing layer whose Mach number is such that the 
corresponding vortex sheet is marginally stable according to Miles' criterion. It is 
shown that the amplitude of the mode obeys a nonlinear integro-differential equation. 
Numerical solutions of this equation show that, when the obliqueness angle is less than 
7 t/ 4, the effect of the nonlinearity is to enhance the growth rate of the instability. The 
solution terminates in a singularity at a finite streamwise location. This result is 
reminiscent of that obtained in the vicinity of the neutral point by other authors in 
several different types of flows. On the other hand, when the obliqueness angle is 
more than jr/4, the streamwise development of the amplitude is characterized by a series 
of modulations. This arises from the fact that the nonlinear term in the amplitude 
equation may be either stabilizing or destabilizing, depending on the value of the 
streamwise coordinate. However, even in this case the amplitude of the disturbance 
increases, though not as rapidly as in the case for which the angle is less than ir/4. 
Quite generally then, the nonlinear interaction between two oblique modes in a 
supersonic mixing layer enhances the growth of the disturbance. 
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1 . Introduction 

Recent interest in the linear and nonlinear stability of supersonic boundary layers 
and free-shear flows has been sparked by the occurrence of these flows in the external 
and internal aerodynamics of high-speed airplanes and their propulsion systems. While 
the subject of the stability of supersonic flows is very intricate, there are (at least) three 
rough rules whose general validity has been confirmed in many independent studies: 
First, the instability of a shear flow becomes more benign with increasing Mach 
number; second, the most unstable modes are oblique; and, finally, there are several 
modes of instability— some of these may appear or disappear for certain ranges of Mach 
number and temperature. Because of this intricate behavior (which can be calculated 
numerically), it is desirable to obtain some theoretical understanding of certain limiting 
cases in order to provide a framework for a physical understanding. 

The present study deals with the (nonlinear) stability of a supersonic mixing layer 
in the long-wave limit. The classic work in this general area is due to Miles (1958) and 
Fejer & Miles (1963). who examined the small-disturbance (i.e. linearized) motion of a 
plane vortex sheet that is subjected to two- and three-dimensional disturbances, 
respectively. They found that the vortex sheet becomes marginally stable above a 
critical supersonic Mach number, and this Mach number is inversely proportional to 
cos 0, where d is the propagation angle of the oblique mode. In addition, there are 
three of these marginally stable modes, which we label according to their phase speeds 
as the fast, slow, and intermediate modes. These classic results of Miles are clearly 
consistent with the previously stated rough guides for supersonic instabilities. 

On the other hand, the numerical studies of Jackson & Grosch (1989) indicate that 
the slow and fast modes are actually unstable in a mixing layer of finite (but small) 
thickness, with the growth rates depending on the square of the small wavenumber in 
the long- wave limit, while the frequency is directly proportional to the wavenumber 
itself. Therefore, the time scale associated with the exponential growth of the 
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disturbance is very much longer than its period of oscillation. Thus, asymptotic 
analysis (e.g. multiple scales) can be used to describe the instability of the base flow. 
Furthermore, the effects of small nonlinearities may be incorporated into our analysis in 
order to affect the growth of the fundamental disturbance through cubic interactions. 

The present study is very different from that of Artola & Majda (1987). They 
examined the resonant response of a supersonic vortex sheet to an incident acoustic 
wave and determined the conditions for the evolution of a kink mode (i.e. a 
discontinuity in the slope of the vortex sheet). The principal difference is that Artola 
& Majda were not concerned with the details of continuous base velocity and 
temperature profiles; their assumed wavelengths are so long that, for all practical 
purposes, the mixing layer appears to have negligible thickness. That is, they were 
interested in the nonlinear evolution of the vortex sheet as if it were a sharp contact 
discontinuity separating the two uniform external streams. This, of course, is not the 
case for the problem here, in which the wavenumbers (or frequencies) are assumed to 
be somewhat higher so that the nonlinear evolution of the disturbance is entirely 
dictated by the details of the continuous base flow profiles, especially in the vicinity of 
the critical level. At this level, the convection speed of the instability mode matches the 
local fluid speed. 

It should be clear from the preceding remarks that the importance of the 
nonlinearities in the present study has its roots in the critical layer where the instability 
mode is approximately steady (in a reference frame moving at the convection speed). 
Thus, the familiar nonlinear terms in the usual substantial derivative cannot be ignored 
in this layer with respect to the unsteadiness; in essence, we have a nonlinear critical 
layer. 

The importance of critical layers in the stability of shear flows has been recognized 
for a long time (Maslowe 1981). Usually critical layers arise for instability modes near 
the neutral wavenumber (or frequency), but they can also arise at low frequencies in 
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incompressible or compressible wall-bounded flows (Goldstein. Durbin & Leib 1987; 
Gajjar 1993). However, they cannot occur in an incompressible mixing layer, where 
the instability is of the Kelvin-Helmholtz type, because both the growth rate and 
frequency are proportional to the (small) wavenumber. 

The asymptotic form of the nonlinearity in the critical layer may manifest itself in 
several canonical forms, depending on the precise balance between the (slight) 
unsteadiness, the nonlinear terms, and (possibly) the viscous terms in the equations of 
motion. One of these canonical forms is the classical Landau-Stuart cubic term in the 
evolution equation for the amplitude of the fundamental (Drazin & Reid. 1981, p. 370). 
Other forms, probably much less familiar but equally important and ubiquitous, are the 
strongly nonlinear phase-jump condition discovered by Benney & Bergeron (1969) and 
subsequently developed by Stewartson (1981), Goldstein et al. (1987), and others, and the 
Hickernell (1984) cubic integral term in the amplitude equation. 

In order to describe the evolution of a symmetric pair of interacting oblique 
instability modes on a shear flow, the integro-differential equation of Hickernell, or its 
variants, appears to be the most prevalent (Goldstein & Choi 1989; Goldstein & Leib 
1989; Goldstein & Lee 1992). In this case, the nonlinearity in the critical layer occurs 
in a hierarchical fashion: The lowest-order interaction of the fundamental mode with 
itself generates the mean flow distortion and the first harmonic, and a subsequent 
interaction of these newly generated disturbances with the fundamental reproduces the 
latter at cubic order in the (small) amplitude. For a properly posed physical problem 
involving spatial instability, each member of this hierarchy of equations is solved with 
vanishing boundary conditions imposed far upstream. The respective solutions are 
expressible in terms of integrals of the yet-unknown amplitude, and the products of 
these integrals, caused by the nonlinearities, ultimately yield the integro-differential 
equation for the amplitude. 
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Since the flow outside the critical layer is unsteady, the nonlinear terras are 
relatively unimportant here, so the flow is governed by linearized equations. Because 
of the very tedious and distracting algebraic complexity that arises in the nonlinear 
problem, we first obtain the linearized solutions. These solutions will be valid 
everywhere (including the critical layer at lowest order). Then, in §7, we calculate the 
nonlinear terms in the critical layer and determine the effects of these terms on the 
amplitude equation. The key result is (44). 

The outline of the paper is as follows: The problem is formulated in §2, and the 

scalings associated with the asymptotic analysis are introduced. It is shown that the 
cross-stream behavior of a low-frequency instability mode is determined by a triple- 
layer structure in the sense of matched asymptotic expansions. The relevant solutions 
in each of these layers are derived in §§3, 4, and 5; this discussion culminates in 
amplitude equation (31^3) for the linearized problem. 

In §6, the amplitude scaling for the nonlinear problem is introduced. This is done 
by requiring that the slight unsteadiness of the fundamental (i.e. its weak instability) be 
balanced by the nonlinear convective derivative, as well as other nonlinearities in the 
equations of motion. The resultant nonlinear equations in the critical layer are solved 
in §7 and in Appendix B. The net result is the appearance of one extra term in (31a); 
the normalized form of this term is given by the integral in (44). This is the final 
theoretical result--it includes cubic-type integrals of the amplitude over the entire life 
of the instability mode. 

In closing this section, there are three points that we wish to emphasize. First, the 
backbone of the nonlinear theory presented herein is the linear theory. The latter 
provides a linear differential equation for the amplitude. Once the nonlinearities in the 
critical layer are incorporated into the analysis via the distinguished scaling of the 
(small) characteristic amplitude with respect to the (small) characteristic frequency, the 
amplitude equation is changed into a nonlinear integro-differential equation. 
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Second, the lowest-order solution in the critical layer may be written as the 
corresponding solution in Goldstein & Lee (1992), denoted by GL, plus an extra term 
that is proportional to a parameter A. In other words, if we (artificially) set A - 0. we 
recover the results of Goldstein & Lee for the flow field of the oblique modes. 
Subsequent products of our solutions always reproduce GL plus some extra terms that 
are proportional to suitable powers of A. In order to simplify the presentation of the 
nonlinear results, we only give these extra terms in Appendix B. Implicit is the 
understanding that the GL terms are available in the literature; in fact, both authors of 
this paper derived all terms independently. Agreement was obtained on all fronts, 
including with GL (Lee. private communication). 

Finally, we establish that a nonlinear integro-differential equation describes the 
evolution of the amplitude at low frequencies. It is well known that a similar equation 
also holds in the vicinity of the neutral point. Since the growth rates are (numerically) 
small at large Mach numbers, it appears plausible that this type of nonlinear equation 
should hold for the entire frequency range. Thus we speculate that the qualitative 
behavior of a symmetric pair of oblique modes may be understood from an asymptotic 
analysis of this type. 

2. Formulation of problem 

Consider a unidirectional, inviscid, and parallel shear flow with velocity U(y)\, 
pressure equal to unity, and density and temperature denoted by pQ(y) and 7 ’q(>’), 
respectively. Arbitrary perturbations in this base flow satisfy 


continuity: 

^ + u • Vp + pgV • u + ^^Po + pV • u - 0 , 

(la) 

momentum: 

^ + U • Vu + ^U'i j-r . 

(1^») 

energy: 

^ + u • V/> + 7(1 + p)V • u - 0 . 

(Ic) 

state: 

p ~ pTq + PqT + pT . 

(1^0 
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where u - (u, o-, w), p, p, and T are the perturbations in velocity, pressure, density, and 
temperature, and the prime denotes differentiation of a function with respect to its 
argument (here, it is y). The effects of dissipation and gravity are ignored, and we use 


D_ 

Dt 



He) 


to denote the time derivative following the base flow. The geometry is defined in 
figure 1 ; x - (x, >», z) is a Cartesian coordinate system, t denotes time, and V - i(3/9x) + 
y(d!dy) + k(3/3z) is the gradient operator. The equation of state for the base flow 
implies PqTo - I. 

All variables are considered nondimensional unless otherwise indicated. Lengths 
and velocities are normalized by a reference length and a reference speed i/ref)* 

respectively. These quantities measure a characteristic thickness and a characteristic 
speed of the mixing layer. For the purpose of developing the analysis, there is no need 
to be more specific about these quantities, whose ratio, provides a transit time 

for the normalization of time, t. Density, pressure, and temperature are normalized by 
their (dimensional) values in the upper stream at - +oo, but the undisturbed pressure 
is actually a constant (set to unity) throughout the mixing layer (as noted earlier). Thus, 
Pq(+oo) = ro(+oo) =» 1, and ro(-oo) - const, actually represents the temperature ratio 
across the mixing layer. We use the notation (*)+ in several different ways: when 

applied to the base flow it denotes the value of a quantity, (•), at y » ±oo; sometimes this 
value may be unity (e.g. Pq^) but. in order to subsequently write our formulae 
economically and "symmetrically'’ with respect to the variables in the upper and lower 
streams, we do not evaluate (*)+ explicitly. 

In view of our nondimensionalization, M - where ^he 

(dimensional) unperturbed speed of sound in the upper stream, is the corresponding 
temperature, and ^ is the (dimensional) gas constant. Clearly, A/ is a characteristic 
Mach number for the flow. Our compressible mixing layer is assumed to be an ideal 
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gas with constant specific heat capacities, whose usual ratio yields the isentropic 
exponent 7 » const. • 1.4. Note that the energy equation (Ic) is actually a combination 
of the first law of thermodynamics, the definition of entropy, and the equations of 
continuity and state. This combination is convenient to use because, in principle, {la-c) 
provide p, u, and p, and (1^ yields T. Therefore, by using this formulation, we obtain 
a certain amount of decoupling among the equations. 

To be specific, let the perturbations in the parallel shear flow be caused by a 
small-amplitude (of characteristic size 0 < e « 1), three-dimensional disturbance. Let 
this excitation occur far upstream (essentially at x - 00 ) in the form of two linear 
instability modes of equal (complex) amplitude and of identical (real) frequency, but of 
opposite spanwise (i.e. z) wavenumbers, which we will assume to be real. Suppose 
further that the frequency of excitation is small. Without loss of generality, we take 
Ai/ - - f/ > 0 and, in order to ensure that the mixing layer is convectively unstable, 

is considered sufficiently positive (Huerre & Monkewitz 1990). 

In order to continue with our formulation, we invoke a number of facts on the 
behavior of low-frequency or long-wave instability modes in a supersonic mixing layer 
at high Mach numbers. An a priori acceptance of these facts by the reader facilitates 
the introduction of the scales for the asymptotic theory, the multi-layer structure of the 
flow in the cross-stream (i.e. y) variable, and a discussion of the cumulative role of the 
nonlinearities. However, the validity of these facts should become clear as we develop 
our analysis; for the time being, we use the corresponding vortex-sheet results to 
motivate the discussion. 

First, an important similarity parameter that determines the growth rate of an 
instability mode in a compressible mixing layer is the Mach number: m - A/ At/ 
(Jackson & Grosch 1989; Balsa & Goldstein 1990). Second, when m is larger than a 
critical value, say the mixing layer is marginally stable in the classical vortex- sheet 
limit (Miles 1958; Fejer & Miles 1963). The actual value of depends on the 
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temperature ratio across the mixing layer and the spanwise wavenumber of the 
disturbance. For example, for a two-dimensional mode and a temperature ratio of 
unity. - 2-\/2. In our analysis, we focus attention on mixing layers with m > 
we call these mixing layers "highly" supersonic, although we are not examining the 
"hypersonic limit" m -* co. Third, these marginally stable vortex-sheet modes convect 
downstream at speed - 0(1) (of course. U_ < < U^) such that \U^ - U^\ is 

supersonic (in the usual aerodynamic meaning of this term) with respect to the 
upper/lower streams, respectively. 

Consider now a low-frequency excitation of a highly supersonic mixing layer. To 
be precise, let 5 - 0(1) be the vorticity thickness of the layer and X be the characteristic 
wavelength of the instability mode: we shall examine the nonlinear spatial evolution of 
a pair of oblique modes when 5 A “ 0(a) « 1, where a - Here cj denotes the 

radian frequency of excitation, and we may think of a as the streamwise wavenumber. 
These oblique modes are actually unstable, with growth rates of 0(a^). Therefore, they 
appear marginally stable only in the classical vortex-sheet limit (i.e. as u -► 0 in some 
sense; see §3). 

These considerations imply that the complex phase speed, c, of a low-frequency 
instability mode is of the form 

c - + ••• • (2) 

where (•) stands for some complex number. Therefore, the thickness of the critical 
layer is 0(a). This layer is centered on the critical level y y^, where U(y^) - U^. We 
emphasize that is given by classical vortex-sheet theory (see (7a)) so that the last 
equality defines y^. Since, in a mixing layer. U(y) is a monotonic function of y, each 
mode has one critical layer. We note that even if this were not the case, it would be 
possible to extend our analysis at the expense of algebraic complexity (see the recent 
work of Wu. Lee & Cowley 1993 on the Stokes layer). Recall that a supersonic mixing 
layer may support two instability modes (the slow and fast modes; see Jackson & 
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Grosch 1989); since these modes have separate critical layers, the nonlinear interaction 
between these modes is too weak to be considered in our analysis. We shall justify this 
remark a posteriori in §6. Thus, the nonlinear interactions that we study arise from 
modes that share a single critical layer. 

In order to describe the structure of the perturbations in the cross-stream 
coordinate, y, the method of matched asymptotic expansions is used. Because of the 
presence of three length scales (X » 5 » a), the various layers in the flow are defined 
as follows: 

y - <9(X) - 0(cr'*) , outer layers (essentially the external streams). 

y - 0{8) - 0(\) , main layers (most of mixing layer where U' - AU/dy + 0), 

y - y,. ~ 0(<j) , critical layer (thin layer centered on critical level). 

Thus, we have a triple-layer structure for the disturbance; the geometry of these layers 
is shown in figure 2. 

The streamwise evolution of the unstable perturbations and the cumulative effect of 
the nonlinearities are described by the method of multiple scales. In view of (2), 
introduce 

I - a{x - U^t) , x^ - a^x , (3a. b) 

for the "fast" and "slow" streamwise variables. Note that on the length scale of the 
mixing layer (i.e. its thickness), even | is a "slow" variable; this is because we have 
long waves, which produce a gently varying disturbance in the x-direction. The 
streamwise oscillations of the instability wave occur on the space variable and the 
slow exponential growth and the cumulative nonlinear interaction between the oblique 
modes occur on the (even) slower variable X| . In keeping with the ideas of multiple 
scales, we note 
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and, for spatial instability, 



+ or^ 






(3c) 


Od, e. /) 


where 0 < erf « 1 is the spanwise wavenumber. The propagation angles of the two 
oblique modes are ±ra/i‘* (crf/cr) - t?(l). Note that 9/3{, 9/9x| . and 9/9f are scaled 
operators that produce an 0{\) effect when applied to the low-frequency perturbations 
in the mixing layer. These considerations are valid in all layers; the correct scaling of 
the cross-stream variable, y, depends on the specific layer under discussion (see 
figure 2). 


3. Outer layers 

The relevant length scale in the outer layers is the wavelength. This implies that 
the scaled cross-stream variable is T = 0(y/\) - ay. The outer layers are defined by 
the condition Y -0(1); in these layers, U - U(y) - U(Y/a) - i/* - const, as cr -► 0, under 
the assumption that U{y) approaches exponentially quickly as ±oo. Similar 

remarks hold for the unperturbed density and temperature profiles. Thus, the outer 
layers are the two external streams in which the base flow is completely uniform; this 
is essentially the vortex sheet limit. The space derivatives d/dx, d/dy, and 9/9r are 
equally important. 

All perturbation quantities (velocity, pressure, etc.) expand as 


velocity: 

u - (u, w) - + oru^'^ + ...), 

(4fl) 

pressure: 

p - ia(p'^^ + + ...), etc.. 

(46) 


where the superscripted quantities are 0(1). Since we are (at the moment) dealing with 
linearized theory (see the discussion at the end of the "Introduction"), the factor ea on 
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the right-hand sides of (4) has no significance; it is introduced here for convenience in 
anticipation of the nonlinear problem, for which we need a characteristic amplitude 
e « 1. This is determined in §6. 

Using (3c. d), we may express our substantial derivative (le) as 

After substituting (4, 5) into (\a-d) and collecting terms with like powers of cr, we find 
that the lowest-order perturbations are given by a slight generalization of the Ackeret 
solution of classical supersonic aerodynamics. Thus, 


- 7 A/ 2 ^(X| )exp[/({ - q^Y + f)] . 


where i - 


and 





7^2 

7b. ‘Zt 


Vt-Uc 

7 A /2 

<7b. 


U^-U, 

7 A /2 


- {fii - P)l/2 > 0 , 


(6fl) 

( 66 ) 

(6c) 

(64) 

(6c) 

( 6 /) 


We recognize /3+ > 0 as the Prandtl-Glauert factor and recall that q* will be real 
above the critical Mach number (Fejer & Miles 1%3), so that the instability wave must 
satisfy a radiation condition at T « ±c». This has been incorporated into the solution 
represented by (6) in order to satisfy causality via "outgoing" waves. 
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We may think of our lowest-order solution in the outer layers as a "marginally 
stable" vortex-sheet mode whose complex amplitude. yCix^), actually changes very 
slightly (by 0{a)) over a streamwise distance of one wavelength. In order to get the 
solution in the upper (or lower) streams, we use the + sign (or - sign) throughout (6). 
Note that the amplitude yC is identical in both streams (i.e. yC^ - yC_ - yC). This 
comes from the asymptotic matching of the solutions through the intervening layers, or 
simply from (6a) and the physical observation that the lowest-order perturbation 
pressure must be continuous "across the vortex sheet." 

In order to make a closer connection with classical supersonic aerodynamics, we 
observe that, at lowest order, the outer field is produced by a "wavy wall" whose 
vertical displacement may be estimated. Using (4a). (5). and (6c), we find the 
displacement 

f )2 )«xp[/({ + f)] • (6g) 

Clearly, this "wall" has a unique displacement; therefore. 


T’o?. 

(V, - i/,)2 " (U. - f4)2 • 


Equation (7a) serves as the dispersion relation whose solutions yield the convection 
speed, t/p, of the instability modes as a function of the spanwise wavenumber and the 
properties of the external streams. Although it is possible to obtain (7a) from the 
asymptotic matching of the solutions across the various layers, it is helpful to have this 
result at hand in order to simplify the subsequent algebra. In this regard, define 


A* - A 


(t4 - ^4)2 


> 0 . 


(lb) 


We call A the vortex-sheet displacement parameter; its effect is felt in the other layers 
as well and shows up in the critical layer in the modification of the solution. 
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As expected, the solution at the next order is considerably more complicated: the 
homogeneous solution is essentially in the form of (6) with 0* playing the role of ^'t, 
and the particular solution has a secularity in Y. This secularity is well known in low- 
frequency theories and is completely unimportant for our purposes. In any case, it may 
be eliminated by further subdividing the outer layers into additional layers; the 
secularity is due to the conversion of a slowly decaying exponential in Y (due to the 
fact that the wave speed, c. has a small imaginary part) to its linear approximation by 
its Taylor series. 

The second-order perturbations account for the fact that the outer flow is slightly 
unsteady in a reference frame that translates with velocity V^\. They are given by 




- 7 A /2 - 5r*r.'t:'(x,)i>expW5 - q,Y + f)] . 


m 


. - 


Ut - u. 


+ u'^^u -^’(JC,)exp[/(| - q^Y + f)] f 


m 


(I) 


a - U, 


(I) 


P Z(iA+f2f4) 


7 A /2 qliU. - U^) 


^'(^l)exp[/(| - q^Y + f)] 


(8c) 


where 




tTr 




a - ^ 


( 1 ) 


P iU 

^ _-<^'(x,)exp[/(^ - q.Y + D] 


m 




(8e) 


and - d-^/dX| . At this stage. ^3+(xi) are two arbitrary functions, which are needed 
to carry out the asymptotic matching. 

From the observation that the flow is homentropic in each of the outer layers 
(albeit with different entropies), the perturbation density and temperature at the two 
leading orders are 
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P 


(j) 


a) 

p 



(9a) 


7 ^) . 2 ^ , ( 9 ^,) 

for y - 0, 1. 

This concludes our discussion of the flow in the outer layers. This flow is 
extremely simple: it is essentially a small-disturbance and slightly unsteady flow. As 
we shall see, the introduction of nonlinearities requires no modifications in the solutions 
at the orders considered. In other words, the disturbances in the outer layers are 
governed by linear dynamics to the required order of accuracy. Thus, superposition 
remains valid and, because of this, we have focused our attention on only one oblique 
mode, whose propagation vector points into the first quadrant of (x, z) space. 
Furthermore, we have written the solutions in complex form, with the usual convention 
in mind that the physical solution is given by the real (or imaginary) part of these 
complex expressions. For the oblique mode whose propagation vector lies in the second 
quadrant, replace t by (-jf) and f by (-f). We now turn our attention to the main layers, 
in which the gradients of the base flow profiles begin to play a role, although nonlinear 
effects are still unimportant. 


4. Main shear layers 

The relevant length scale in the two main layers is the characteristic thickness of 
the mixing layer. This implies that the scaled cross-stream variable is y/6, or simply y, 
since 5 - 0(1). The main layers are defined by the condition y - 0(1), although the 
critical layer, sandwiched between them, is excluded. We must account for the 
dependence of U, ^0 y these layers. 

The perturbation quantities expand differently than in the outer layers. This 
signifies completely different physics. As we have seen from the outer solutions. 
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disturbances of 0(ea) produce a vortex-sheet displacement of 0{ea/a) - 0{e) » 0(ea). 
Such a displacement of a material volume in a non-uniform base flow causes 
perturbations of 0(e). Therefore, 

u - e{iP^ + + ...) , (lOa) 

with similar expansions for the perturbation density and temperature. On the other 
hand, since 8fK - 0(a) « 1 . the main layers are "thin" so that the pressure is 
impressed on them by the outer layers. Thus, 

p - + ...) (10i>) 

and 

V - + ...) (10c) 


since, in a thin layer, - 0{u SfK). Finally. 

,.s(0) .^(1) , 

w - fff(w + ffW + ...) 


( 10^0 


since spanwise velocity perturbations must originate in spanwise pressure gradients (the 
previously alluded to displacement effect is absent since the base velocity is along the 
x-axis). The caret indicates that a variable is written for the main layers, and the 
superscript designates the place of a term in the expansion. As usual. ( )(°) and ( )0) are 
0(1). The ellipses in (10) stand for higher-order terms that are not needed in this 
study. Using (3c. d), we express our substantial derivative, (le), as 


_D 


dt dx 


- o< 


{U - U,) 




( 11 ) 


and. after substituting (10) and (11) into the linearized form of (1). we may obtain a set 
of equations for ( )(01 and (•)^*^- 

The lowest-order equations are simple to solve. For example, dp^^/dy - 0 so that 
^(0) is a function of |, f, and X| . This function is determined from matching with the 
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pressure in the outer layers. Thus, 

yC{xi)exp[i{i + f)] . (12a) 

and, from the span wise momentum equation, we get 

IT’ 

Note that (124) matches with (6</); it also blows up as >> ■* This establishes the need 
for the critical layer and a rescaling of the solutions there. This is discussed in §5. 
Similarly, by solving the remaining equations of motion, we obtain 


(f f, X, ) , 

(12c) 

- - rbcfic. f. X,) . 

(\2d) 

p*°^ - - Po<7f\f f. X,) . 

(\2e) 

- (U - f, X,) , 

(12/) 


where (7^°^ are two arbitrary functions of the indicated arguments (i.e. independent of 
y). The ± signs refer to the values of these functions in the two main layers, one above 
and one below the critical layer, and the subscript | on denotes partial 

differentiation with respect to Also, U' = AUjAy, T'q - dT^ldy. and pg = dpg/dji^. 

The asymptotic matching of (12/) and (6c) immediately establishes 

- <7^°^ - -/A^(x, )exp[;({ + f)] . (13) 

where A is defined by (lb). Note that (12c, d, e) are perfectly consistent with our 
opening remarks of this section; Eulerian perturbations in the u-velocity. temperature, 
and density are proportional to the corresponding base flow gradients and arise from 
the displacement of a material volume in the nonuniform base flow. 
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We are also interested in the solutions at the next order; it is precisely here that 
the main-layer expansion further breaks down as >» -*• ^-momentum equation 

may be solved to yield the pressure at second order 


where 


'(1) 


- 


„(0) 


ry 


y + 


[r*0;) - l]dy 


±00 




rttv) 


U(y) - Uc 
44 - 44 

T’oO^) 

^0, 


(14a) 


(144>) 


and are two arbitrary functions of their indicated arguments; they represent the 
homogeneous solution to the boundary layer equation d/dy » 0. It is possible to 
determine them by matching the pressures in the main layers to those in the outer 
layers. We find 

pP - 7A/2 0j(x,)exp[/(f + ?)] . (14c) 


where was introduced in (8a) as part of the homogeneous solution in the outer 
layers at second order. 

The span wise gradients of the pressure given by (14a) drive the spanwise velocity 
component. The corresponding equation of motion may be solved for 

>^(0 

■ - 57^ -CU,)exp[/« . OJ - ^ ^ . (15) 

Note that the first term on the right-hand side of (15) arises from the variation of wW 
on the slow scale X|, whereas the second term is due to spanwise pressure gradients. 
The spanwise velocity components in the outer and main layers automatically match. 
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At second order, the streamwise and cross-stream velocity components are most 
easily obtained from the energy (Ic) and x-momentum (li») equations. The algebra is 
quite involved; therefore, we provide an intermediate step that documents the coupled 
equations for and 


^ dy 3X] 7 


iU-U,) - f/ _ - ^ . 


(16a) 




After eliminating 9u<*)/9| from (169) via (I6a), substituting the lowest-order solutions 
(12) into the forcing terms, and solving, we obtain 


-d) 




+ (U - f. ;c, ) , 


(17a) 


where 


F.M, y. r. X,) - J 


ry r 


(1 + |2) 


±cx> 


rtfv) 


1 


dy - yQ± -^(jci )exp[/({ + f)] 


+ 9i^'(xi)exp[/(| + f)] , (179) 

and (j'P is an arbitrary function of its arguments; the last term in (17a) represents the 
homogeneous solution. 

The streamwise velocity component follows from (16a) or (169) after eliminating i^(*) 
via (17a). We find 

^ y' f’ “ u^-U + f)] 


- f. X,) . 


(17c) 
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The asymptotic matching of the cross-stream velocity component to that in the outer 
layers yields 

- -/AS3i(x,)exp[/(f + ?)] 


+ 






+ 



\yC'{xi ) exp[/({ + f)] . 


( 17^0 


As mentioned before, the solution in the main layer breaks down in the vicinity of 
the critical level where {U - U^) is small; see for example {\2b). In order to repair the 
resultant disorder in the expansion, the critical layer is introduced. The behavior of the 
main solution at the edge of the critical layer is discussed more fully in Appendix A; 
such a discussion is essential in order to match the solutions in the critical and main 
layers. 


S. The critical layer 

The critical layer is centered on the critical level, y ~ and its thickness is of 
0(a). Therefore, introduce the critical layer coordinate 


with T) - 0(1), and expand the base flow in its Taylor series about the point y 
We find, using (le). 


(18) 


>'c- 


D 3 r; 3 2 

d7- 


+ 0(a2)|. 


where .2^®) and are two linear operators defined by 




31 • 


^■) 


,/■ _9_ ^ „ a. 

3x, 2 ^ di 


(19fl) 


( 196 ) 

(19c) 


The notation ( )g means that (•) is evaluated at y ~ y^. 
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Based on the discussion presented in §4, we find that the perturbation quantities 
expand as 


MO) A(l) s 

u - €(u + au +...), 

(20a) 

MO) M^) ^ 

w - e(w + (xw +...)» 

(206) 


with similar expressions for the density and temperature, whereas 

p - ea(p^^^ + + ...) (20c) 

and 

l>- m . (20d) 


Note that (20c, d) are entirely consistent with the fact that the critical layer is a thin 
boundary layer (in a generalized sense of this term), and the form of (20^>) expresses the 
fact that the spanwise velocity component in the critical layer is 0{\/a) larger than that 
in the main layers (see 12(6)). 

The double caret indicates that a dependent variable is written for the critical 
layer, and the superscript denotes the place of a term in the expansion. As usual, ( )(°) 
and ( )0) are 0{\), and they depend on the critical layer coordinate t) such that 3/3 tj of 
these terms are also of 0{\). 

We now use (18)-(20) in the equations of motion, (1), in order to derive a sequence 
of equations for the flow in the critical layer. The simplest of these comes from the y- 
momentum equation and yields 


dr\ ~ dif 


( 21 ) 


so that pb) . /<J)({, f, X|) U - 0, 1); the actual function form may be obtained from 
matching with the main layers. Thus, using (12a) and (14a. c). we find 


- yM^Mxi)exp[i(( + f)] . 


(22a) 
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Vc + 


[rtCv) - l]d;; kY'(jC|)exp[/(| + f)] 


±00 


+ 7/\/2£gj(x,)exp[/(C + f)] . 


{22b) 


Since there must exist a unique p(') that is independent of the label (±), the 
solvability of {22b) implies 


where 


^+(^ 1 ) - ‘Q^~^{xy) - S.{x^) - ^Q..'C{x^) , 


Q, 


const. 


Qi 


rye 

[riCv) - i]d>^ ► . 

J ±00 


(23a) 


(23b) 


Now recall from (8) that <55^ are associated with the homogeneous solutions in the outer 
layers at (9(cr). These solutions merely reproduce the lowest-order solutions; in order to 
assign a unique gauge to the amplitude we set * 0. Thus, we interpret as 

the amplitude of the instability mode (to the required order of accuracy) in the slow 
external stream, and (23a) defines . 

TTie lowest-order spanwise velocity component is essentially generated by the 
f-derivative of ; we solve the leading-order z-momentum equation to obtain 

JT^ ^q{V. jc,)exp[/(| + f)] , (24a) 

where 

^i) - ^i)exp<-/0(T;x,) , {24b) 


^i) 


r^l 

(J| - Xi)"^(X|)eXp(/Q£TJX|)dX| , 

J_oo 


/I « 0. 1, 2 


(24c) 


and Q! = U\!U^. This apparently cumbersome notation involving If'n and /„ turns out 
to be extremely useful for representing the hierarchy of solutions of the linear and 
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nonlinear problems. Therefore, we introduce it at this stage of the analysis. Note that 
the asymptotic behavior of {n - 0, 1, ...) for large \rj\ may be obtained by 
successive integrations by parts. For example, 

^i) - (25) 

^ ‘ lOLr] {iOLr\Y 


under the (correct) assumption that vanishes exponentially as X| -oo. Similar 

expressions may be obtained for ^ 1) whose leading term is 

The x-momentum and energy equations for 3^®) and are coupled. However, 
an advantageous decoupling is effected if we first solve for du^^'^ldr). The relevant 
equation (essentially for z-vorticity) is 


9xi 


a 


dfj 


Oil 




(26a) 


which expresses the fact that a (z - z) normal strain stretches a spanwise vortex line 
and thereby changes the corresponding vorticity. The solution of (26a) may be 
integrated immediately with respect to rj to yield 


“ - W^o(n. ^i)exp[/($ + f)] + /Af/;.-.^(x,)exp[/(| + f)] , (266) 

where the last term in (266), the "constant" of integration, arises from matching (266) 
and (12c). Finally, having determined and S(®), we obtain from the lowest-order 
x-momentum equation 




iTo (1 + f2) 

— -/^(x,)exp[/({ + f)] - i^.U^\.yC'{x^) 


+ /QiT>^(x,)]exp[/({ + f)] • 


(26c) 


Note that (26c) automatically matches with (12/). 

At this point in the analysis, v//(X[) is still an arbitrary function of its argument; in 
order to determine the dependence of the amplitude on the slow variable X|, we must 
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calculate the 0(a) corrections to the critical layer solution. These are denoted by 
in (20); specifically, we need 30 ) in order to derive the amplitude equation. 

The spanwise velocity component, wO), obeys the inhomogeneous equation 




- - 


i d_ 
tA/ 2 df 




(27a) 


whose right-hand side is known from the lowest-order critical layer solutions, {22a. b) 
and (24a). In order to solve (27a), we observe the identity 



+ iOLI) 


lT,(i?, X,) - - W'o(ij. X,) . 


(276) 


and write the solution 


7j 2 fy.. 1 IT- Q 

w(»expH(? + f)] - ^ u — 

+ ^ (aTg^ - 7’o^)l^o('J* ^i) • (27c) 

It is possible to show that the two-term critical-layer solution for the spanwise 
velocity component matches with that in the main layers in the sense of the asymptotic 
matching principle of Van Dyke (i.e. 2-main of (2-critical) - 2-critical of (2-main)). 
For this purpose, the behavior of the solutions in the main layers near the edges of the 
critical layer is derived in Appendix A. 

As before, it is simpler to solve for the 0{a) correction to the spanwise 

vorticity, rather than directly for The governing equation is 


- f + f/’; ^ - 1/-; 


whose right-hand side contains the vortex stretching and displacement effects, as well as 
that of the baroclinic torque (the last term). 



- 25 - 


If we substitute the previously determined solutions into the right-hand side of 
(28a), we obtain, after a modest amount of algebra. 


exp[-i($ + ^ 




- 2a^ 


^i) - 


ifiotTQ Q_ + Ifia^TQ r} 


+ i/’; ToJWQir,, xO + iT. 


U, 


U",{\ + fi) _ 


U' 


>yC(xx 


+ /At/'c exp(-/aT?X|) ^ — [^(X]) exp(/ai;X|)] 

C 7 X I 


(286) 


whose solution is 


9tj 


ifialQ U" ,2 

exp[-/(? + ?)] ^2 ^l) - jj2 (^c^bc + ^l) 


'To. [(/" (1 + P) T’o 
t/. 1 U' To 


^i) JJ 


fia^To rt 


dx. 


Q. 

+ 77 ^ ^i(ri, X,) + /Af/-; Mxi) . (28c) 


In deriving (28c), we have made use of an identity, which is a slight generalization of 
(27Z>), involving the higher-order 's. 

In order to carry out the asymptotic matching with the main layers, we first 
reconstruct the streamwise velocity component 


U^'^TJ) - S^'^O) 


r ^ 


du 


,?.0) 


dt) 


dt? , 


(29a) 


where, for emphasis and brevity of notation, we only show the tj-argument of 30 (of 
course, the other arguments are |, f, and ). We cannot let tf •* ±oo in (29a) in order to 
determine the behavior of at the edges of the critical layer for the purpose of 
matching, since the above integral is divergent at the upper limit. 
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In order to get around this difficulty, introduce 


Z(v) - 


3t) 


I 'T B 

iAU'l. - -jp— (1 + |2) 'yC(xi) exp[/(? + f)] 


( 296 ) 


where 


B 


Hl. ^ 


U' 


(29c) 


Jy - Xc 


Observe that it is possible to show from (28c) and the asymptotic behavior of that 
Z(t>) is analytic for all (real) t? and is integrable at ±oo. Therefore, we rewrite (29a) as 




a^‘^(0) 


r 

Z(ij)dt? + 

Jo 


fj? 

Z(l?)dTJ 

j ±00 





T B 1 

(1 + fi) logd + jj 2)1/2 |•^(J£;,) exp[/({ + f)] 


(29d) 


and note that the second integral in (29d) is of 0(i?-') for large t/. 

We now perform the asymptotic matching principle of Van Dyke in the form "2 
main of (2 critical) - 2 critical of (2 main)" for the streamwise velocity component. 
With the help of Appendix A. (29d), and (266), we find 


r -yc 

QiVc - ~~ J*(y)<^y 

J±oo 

X yt{x^) exp[/(| + f)] + /?(!. ?. JC|) . (30) 

where /? is a nonessential remainder that is actually known from matching, but whose 
precise form is not needed for the determination of the amplitude equation. There is 
no (±) label on R, Recall that is defined by {\ld) (0± in this equation may be 
obtained from (23a) in which we set 0^ » 0; see earlier discussion) and Jt(y) is 
introduced in (A 4). 


^±oo 


w^'\o) 


Z(Tj)dTj - + AU’ 
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We now have two equations for the single unknown D(*)(0). This is because we 
may match the critical layer solution to the (±) main layers. The existence of S(’)(0) 
(i.e. solvability) requires 


exp[-/(| + f)] 


Z(Tj)dfj - + i;2^(X|) , 


(31a) 


where 


- At/ci/U 


U - fA 


1 - 


1 + fiUJU. 


1 


i/. -f/c 


. 1 + 


(316) 


^2 = mi 


^1+12 

r 

J.{y)dy + 

Q- 

J 


r +00 


(y)dy 



r>^c 

+O0 


[r.(>') - l]d>/ + 

[r,(>') - l]dy> 


~oo 



(31c) 


and r+fv) is defined in (146). 

We interpret (31a) as the evolution (amplitude) equation for For a given 

mixing layer, the constants and i >2 may be calculated, and the integral of Z may be 
obtained (most simply) as follows. 

Since the desired (improper) integral of Z exists in the Riemann sense, it also exists 
in the sense of the Cauchy principal value, lim It immaterial, in 

b“*oo 

principle, which of these two procedures is used to calculate the value of this integral. 
However, from practical considerations, it is simpler to calculate the Cauchy principal 
value because the second term in the braces of (296), being an odd function of rj, 
integrates to zero. Thus, 


fOO pOO pOOp^^^ 

Z(Tj)di 7 - i Z(i})d7j - i - - ihU"^yC{x^) exp[/(^ + f)] [dj) . 

3 -oo J -oo J -O0 L J 


(32a) 
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After substituting (28c) into (32a) and using the identity 


r 

J-oo 


2ffkt 

(-ia)^ I aj 


Mxi ) 


fO 

6(x)djc , 

J-oo 


A - 0. 1. 2 


(323) 


we find 

1 ni To B(l + f2) 

J ool^ ’ ” ■ R ~~^c 

X Mxi) exp[/(| + f)] , (32c) 

where we have assigned the (correct) value. 1/2, to the integral of the delta function in 
(32^). Note that (32i>) is a purely formal result that may be established by the use of 
generalized functions (Lighthill 1978); the two possible values for the integral of the 
delta function are 0 and 1 (whenever the upper limit takes on values 0“ and 0*. 
respectively). 

Although it would be tempting to conclude from these remarks that the value of the 
integral on the right-hand side of (32b) is 1/2, such a conclusion can only be arrived at 
by a more careful procedure. This involves the ideas of "phase jump" associated with 
the logarithmic branches that arise in the stability of shear flows with critical layers 
(Maslowe 1981). As is well known, this phase jump depends on the algebraic sign of a 
» hence the presence of | a| in (32c). Our final remark deals with a 

generalization of (32b) that is also needed: When the exponent of is less than the 
value of k on the subscript of W^, the integral vanishes. Therefore, only the first three 
terms on the right-hand side of (28c) yield a nonzero contribution. 

If we now substitute (32c) into (31a), after invoking (32a), we indeed arrive at an 
amplitude equation for ^(X| ). 

In the nonlinear theory, which is developed in the following sections, the amplitude 
equation (31a) suffers a modification. Briefly, the right-hand side remains unaltered, 
but the left-hand side will contain a term additional to that given in (32c). 
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6. Amplitude scaling for the nonlinear problem 

In our discussion so far, the characteristic amplitude e « I has not played a role 
because the problem has been linearized. We now establish the distinguished scaling, 
e = £(a), which will enable us to include the effects of nonlinearities in the phase jump 
condition, (32c). 

Of course, as is seen from (1), there are many nonlinear terms in the equations of 
motion. In order to motivate our discussion and the derivation of the scaling, we shall 
use the convective nonlinearity, u • Vu. as a typical term; specifically, we write the x- 
component of (li). 


Dr 


+ u 



- 0 , 


{33a) 


and focus attention on the first two terms of (33a). 

The nonlinear terms first show up in the critical layer because the linear term, 
Du/Dt, is relatively small here. After using (19a) and (3c) in (33a), we find 

+ ...|u + au ^ + ... - 0 . (336) 

Since, from (20a), u - 0{e). the nonlinear "forcing" term in (336) is 0{cre^) and the 
resultant first and zeroth harmonics are of 0{a€^/a^) - 0{e^/a). Note that this is 
significantly larger than what would be obtained by simply squaring the fundamental 
(i.e. <9(f2)). A further nonlinear interaction of these harmonics with the fundamental 
of 0(e) produces nonlinear terms of 0{a(€^/a) - It is at this level in the 

perturbation analysis where the fundamental reappears again for the first time due to 
nonlinearities, and, in order to affect the phase jump condition (hence the evolution 
equation for yC{x^)), we require that the slow aging of the fundamental (on the scale Xj) 
also be present at this level. Thus, from (336). 


ff^ore . 0{e^) , 


(33c) 
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or e - C>(a3/2) por convenience, we set c ■ a^/2 without loss of generality. This is our 
distinguished scaling. 

It should be noted that this type of scaling argument was used by Hickernell (1984) 
in a geophysical problem and by Goldstein and his colleagues in a more closely related 
set of problems involving oblique modes near the neutral point (Goldstein & Choi 1989, 
Goldstein & Leib 1989) and resonant triad interactions (Goldstein & Lee 1992). 

It is somewhat difficult to make a direct comparison between our distinguished 
scaling and those of others because of the more elaborate triple-layer structure of the 
present problem (i.e. outer, main, and critical layers) as compared with the usual 
double-layer structure (i.e. outer and critical layers) found by previous investigators 
for their problems. Nevertheless, the following remarks may be helpful: let eg « 1 be 
the characteristic magnitude of the instability mode in the outer layers. We see from 
(4) that €q - 0{e(x) - 0 (cr^/ 2 ) because of our distinguished scaling. Thus, tr - 
when our result is interpreted in this way, it is reminiscent of the scaling used by 
Hickernell (1984) and Goldstein & Leib (1989). 

There remains to assess the importance of nonlinearities in the main and outer 
layers. By using (11) and (10a) in (33a), we find that the nonlinearities do not enter in 
the main layers to the required order of accuracy. Similar remarks hold for the outer 
layers. Therefore, with the exception of the critical layer, the other layers are governed 
by linear dynamics. Thus, the relevant solutions for the fundamental disturbance are 
essentially available in §3 and §4. 

Based on the above remarks and (20). we find that the perturbation quantities in 
the critical layer expand as 

u - a3/2(3^°^ + + ...) , (34a) 


with similar expressions for w. p, and T, whereas 
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p = ff5/2(^(0) + ^ ^ ) (34^) 

and 

^ - a7/2(^<°) + ctI/ 2^<'/2) ^ ^^1) ^ j _ (34^) 

where the ( )(^/2) terms are inserted in order to accommodate the zeroth and first 
harmonics (which are, for example, in the w-equation of 0(t^ fu) = 0(a^)). The nonlinear 
part of the complete substantial derivative, to the accuracy required herein, becomes 

U • V - (t 5/2 + 0'i/2//'^^^ + 0(ff)j . (35a) 

where 

A + ^^) |- + 1^^) §r , / - 0, 1/2 . (35i) 

Ok dv 

After substituting (34) and (35) into (1) and expanding the base flow in Taylor 
series about the critical point, we obtain a hierarchy of equations for (’)^K j - 0, 1/2, 1. 
The lowest-order solutions, ( )^^^ are essentially given by linearized theory; these were 
already obtained in §5 using complex representation. For the nonlinear problem, we 
consider two oblique modes, with spanwise wavenumbers of (±1) and exponential factors 
of exp[/({ ± f)], respectively. The complex amplitudes of these modes remain ^(X|), 
although it is possible to generalize our analysis to unequal amplitudes. 

After superposing the oblique modes and using the imaginary part of this complex 
quantity to represent the physically real solution, we find from {22a), (24<z), (266), and 
(26c), 

a 7 A /2 Re[/V^ exp(/^)] cos f , (36^2) 


tr 

- -T^ Re[/if'o exp(/'5)] sin f , 
fiTn 

jjS. Re[^f'o exp(/0] cos f - Re[_^ exp(/|)] cos f . 


(36^>) 


(36c) 
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MO) 

^ Re[^ exp(/J)] cos f 

+ Ai/(. Re[^' exp(/|) + iar) X exp(/^)] cos f . (2>bd) 

where Re(-) denotes the real part of a complex number, X - yC(x^), X' - d^(X|)/dX|. 
and Wq - Wq(,t). x,). In order to simplify the notation a bit. we make use of the 
convention that, whenever the argument of a function on the right-hand side of an 
equation is not shown explicitly, it is understood to be the usual field variables x, or 
(i?, X|). etc. For example, x,) will be written as /„ (n - 0. 1, ...). 

It should be noted that, if we artificially set the vortex-sheet displacement 
parameter. A. to zero in (36). we recover (albeit in the current notation) the critical layer 
solutions used by others (e.g. Goldstein & Lee 1992). This comparison points to the 
considerable algebraic complexity that we will encounter in the treatment of the 
nonlinear problem. In order to keep this complexity at a minimum, we shall focus on 
the new terms that come about because A j* 0 (see Appendix B). We might point out 
that, because of the presence of these terms, the cross-stream velocity component, is 
no longer a "constant" in the critical layer but varies linearly with t) (see {36d)). 


7. Nonlinear solutions in the critical layer 

Although the linearized theory suggests that we should solve for and 

integrate this quantity from tj « -oo to tj - +oo in order to find the logarithmic phase 
jump (see (32c)). the analysis of Goldstein & Lee (1992) clearly indicates that a 
significant economy in algebraic manipulations is brought about in the nonlinear case by 
introducing the new variable 


Q 


do dif 


(37a) 


whose fundamental Fourier component. 0, , , is defined by 
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and obeys 


0-0,1 exp[/($ + f)] + ... . 




d 

dxi 


+ lar^ 




dri 


+ (linear terms), , . 


(376) 


(37c) 


The nonlinear term. J/', in (37c) arises from the convective nonlinearities in the 
equations of motion. It is defined by 


^ - yo)(3('/2) ^ ^ ^«/2)(S(0) ^ , 


Old) 


with denoting its fundamental component, 

.r - , exp[i({ + f)] + ... . (37c) 

where the ellipses in (376. e) stand for the other Fourier components. Recall (356) for 
the definitions of and M'/2). 

We note that the derivation of (37c) is straightforward by taking linear 
combinations of the equations for 9S(*)/5ij and dw^^^ldrf, observing that dpO/dri » 0 
(J - 0, 1/2. 1). and invoking the pressure form of the equation of continuity (i.e. the 
energy equation) to eliminate (SDO/S? + 9^(')/9i?) in favor of [-(I9w(‘)/9f + 9w(°V9^i)]- 

The phase jump for the nonlinear problem is obtained most easily by using 
0, ,exp[/({ + f)] in place of du^^^/dr) in (32c). This is permissible since the fundamental 
component of 2)0 has no jump across the critical layer; a fact guaranteed by the 
behavior of the main-layer solutions at the edges of the critical layer (see (A 7)). 

Since 0, , obeys a linear, inhomogeneous equation, the phase jump associated with 
0, , is the superposition of the effects of the two forcing terms in (37c). Insofar as the 
"linear terms" in this equation can be easily obtained from (27a) and (28a), and the 
phase jump associated with these terms has been calculated already, we shall not say 
more about these terms. Our principal focus is to calculate the effect of the nonlinear 
forcing term -'k', , . 
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Before we go any further with some general results, it is advantageous to introduce 
the complex Fourier-series representation of a perturbation variable. 


^n!m (’!• exp[/(n{ + wf)]. j - 0. 1/2. (38) 

m n 


where / may stand for u, v, w, etc. Of course, some of the F's may vanish. Note that 
for a function the Fourier coefficients are denoted by Fn.iii- example, if y » 0 
and / we find from (38) and (36rf) 


-(0) 
^ 1,1 

-( 0 ) 


-( 0 ) 

1,-1 

-( 0 ) 

*^-1.1 


& , 
& . 


(39a) 

(39i) 


where 


^ 7'o (l+P) Ai4, 

& - —^ 4^71 + lariA) 


(39c) 


~ A(0) 

and (•) denotes the complex conjugate of (•). All other ^ vanish. Similar 
expressions hold for and ^n.m- which are derivable from (36c. b) and (38). 

In order to calculate the nonlinear forcing term. 3(J^i i)/9t>. we carry out the 
operations indicated in (37c/) by using the Fourier-series representation of the variables. 
The resultant expression may be simplified slightly by employing the equations 


5(1/2) Mi/2) 

ITo.o - ^ 


^ 2.2 + ^ *^ 2.2 * ^ 


(40a) 

(406) 


and the expressions for the Fourier components of (-)^®^. as obtained from (36). Note 
that the validity of (40a. b) is discussed in Appendix B. After a modest amount of 
algebra, we find 


3-^i.i 

3t) 



k~\ 


(41a) 
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where 


T,~& 


32 i/ 


( 1 / 2 ) 

0,0 


9tj2 


r , - (? 




dri^ 




T. — 

'» 2U. 3t? 


A(i/2) aw^ 

^^ 2.0 dr, 


^ A 

i/c 3 t , 


'i/o.2 ^'O + ^ ^'0.2 


/a:/: 


^ -/f: 


3(/ 


r,('/2) 


2,0 


3tj 


#»(l/2) 
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(413. c) 
(41rf) 

(41e) 

(41/) 

(41g) 

(413) 


and Wq is defined by (243). This particular grouping of the terms is helpful in 
reducing the subsequent algebra (Lee 1992. private communication); they represent the 
relevant "cubic" nonlinear interactions that produce the fundamental. 

Unfortunately the algebraic manipulations do become extremely tedious beyond this 
point— we shall be content with presenting a few of the key steps in Appendix B. 
However, the details of these manipulations may be found in the dissertation by 
Gartside (1995). 

A 

The overall plan is to solve for the (•)C/2) variables in the critical layer (see (34)). 
These are forced by the quadratic interactions of the fundamental with itself; the 
solution provides the Fourier coefficients (e.g. etc.) needed in (413-3). We then 

solve (37c) to obtain the seven contributions to Q, ,. say arising from each of the 

nonlinear forcing terms (3 - 1 7). The nonlinear correction to the phase jump 

condition. (32c). may be found according to our previous discussion in this section. For 
the purposes of illustration, some of the algebraic details are carried out in Appendix C 
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for the simplest of the terms, namely, T-j. The other terms may be handled in a similar 
fashion, as shown by Gartside (1995). 

After a considerable amount of algebra, we find that the total nonlinear correction 
to the phase-jump is 



i.i(t). ;c,)d7j 


ijr T\ sin20 cos 20 


16f/g cos^0 


sgn(f/;.) I yC(xj)dx, [ z|z(Z + ^) 

j -OO J -oo I 


^(xi)Sc(x 


1 + X, - jc,)djc, 


and 


tan 0 = t . 


Z - X| - X!| , 


^ - 5, - X, . 


{42a) 


{42b) 

(42c) 

{42d) 


It is assumed that > 0 in these equations. We multiply (42a) by -4/ exp[/(^ + f)] and 
add this to the right-hand side of (32c) in order to obtain the total phase jump (due to 
linear and nonlinear effects) to be used in the amplitude equation (31a). We shall give 
the canonical form of this equation and the corresponding initial conditions (as X| -* -oo) 
in the next section, after introducing a series of transformations in order to reduce the 
number of parameters. 


8. Discussion and conclusions 

In order to present the simplest evolution equation for the amplitude, we introduce 
scale transformations for the independent and dependent variables: 

^ < 0 - (43a) 


^yC{x^) m r>4(x). 


(436) 
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where t^i and U 2 are defined by (31i», c) and F is a complex constant. Recall that and 
O 2 are determined from the base flow profiles and the underlying linear instability 
theory that defines the critical point. The magnitude of F (i.e. | F|) is chosen such that 
the absolute value of the coefficient in front of the integral on the right-hand side of 
(42a) becomes unity and the phase of F, together with a shift in the origin of the slow 
variable x, is used to eliminate a factor of proportionality from the upstream boundary 
condition. Similar transformations have been used by Goldstein and his colleagues (e.g,, 
Goldstein & Lee 1992); the algebraic details are quite straightforward and will not be 
repeated here. 

Using these transformations, as well as (32c) and (42a, b), we find that the 
amplitude equation becomes 


- iKA - sgn(i' 2 ) sgn(i/g) sgn(cos 26) 


pX, 

/4(x,)dX| 

J-oo 


Z^{Z + ^) + 2^(^ - .^sin2fl^/l(5|)/l(5| + X| - X|)dx| , (44) 


where we have renamed our new slow variable (i.e. (43a) and a shift in the origin) to 
be in order to make use of convenient definitions (42c, d) and to keep the notation 
simple. However, the presence of A (rather than in (44) serves as a reminder that 
this equation is written in terms of the "scaled" variables. The upstream condition 
associated with (44) is 


where 


e-Kx, 


as X 


ff sgn(f/e) 
^2 cos^6 




1 - 


ToU", 


(45a) 


(456) 


This condition expresses the fact that very far upstream the base flow is perturbed by a 
linear instability mode. 
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The solution to the amplitude equation. (44), formally depends on three parameters: 
Re(/0, the coefficient in front of the double integral, and d = tan** f. We shall study 
(numerically) the behavior of /4(X|) as a function of these parameters; some 
representative results are presented below. 

In figures 3 and 4 we show representative results for R^K) - as a function of 
Mach number, m - AUM. for three different temperature ratios. Tq /Tq^ - Tq . The 
velocity profile in the mixing layer is given by 

^ tanh(yH) (46a) 

where 

I . AU ~ - U, > 0 {46b) 

and the temperature profile is specified by Crocco's relation 

1 - T 

^oCvh) “ 1 2^ ' tanhCv'H)] + rrfi[l - tanh2(j;„)] (46c) 

where is the Howarth-Dorodnitsyn variable. It is related to the "physical" cross- 
stream coordinate, y, by the transformation 


y 


•yu 

To(i?)di? 

Jo 


( 46^0 


These profiles were also used by Jackson & Grosch (1989) in their study of the linear 
stability of compressible mixing layers. 

It is seen from figures 3 and 4 that, generally, is of order unity and varies 
roughly linearly with m, for m sufficiently large. The fast and slow modes are defined 
by the inequalities and (/g < i/^, respectively, where is the convection 

speed of the instability mode as obtained from the vortex-sheet dispersion relation, (7a). 
Note that when the temperature ratio is unity, the two curves for are mirror images 
of each other about the horizontal coordinate axis. This is because {U - U^) is an odd 
function of y - • 
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Interestingly enough, the solution of the nonlinear amplitude equation, (44), does not 
depend on in an essential way. To see this, consider the trial solution 


A - . (47) 

which assumes that the phase of the complex amplitude is linear in , After 
substituting (47) into (44) and simplifying, we find 


- I - sgn(i; 2 ) sgn(i/;,) sgn(cos 20) | A(Z + x^)| dZ 

^ J_oo 


mtt)\A{l +x^)\\A{t +Z + x,)|d2 


(48a) 


where the kernel, 5f^(Z, is defined by 

5T(1 - Z{Z(Z + ^) + 2^ (^ - ^in20} . (48*) 

This nonlinear equation for \A\ is independent of ATr. 

Furthermore, for our mixing layer, U'^ > 0 and, for profiles given by equation (46), 
U 2 turns out to be negative (as obtained from numerical evaluations of (31c)). Thus the 
coefficient in front of the double integral in (48a) is +1 according to 6 > 7 t/ 4 . The 
algebraic kernel in (48) is negative over the domain of integration. 

Thus, for Q < 7r/4, the amplitude of the oblique mode monotonically increases with 
coordinate x, . This observation is confirmed in figure 5; the corresponding phase of A 
is shown in figure 6. The latter figure establishes the validity of the trial solution, (47), 
in this particular case. 

The numerical result in figure 5 suggests that |/4(X|)| develops a singularity at 
some finite value of x, . say This observation is consistent with the findings of 
other investigators (Goldstein & Choi 1989, Goldstein & Leib 1989); indeed, near the 
singularity, the dominant terms in (48a) can be balanced by a solution of the form 
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Mxi )j 



(49) 


as X, -* from below. 

In order to establish the presence of a singularity when 9 < ff/4, consider the 
inequality 


fO 

I A{Z + Xi)| dZ 

J -OO 


f z 

55'd I +x,)| I A{^ +Z + x,)| d^ 

J-OO 


^ - 


rO 


dZ 

J -oo 



55"(Z, ^)\ A(^ + Z 


+ X[)| ^d^ , 


(50) 


which holds because |/4| is a monotonically increasing function, and Z ^ 0 and ^ ^ 0 
over the domain of integration. Both sides of (50) are positive. The integral on the 
right-hand side of (50) may be reduced to a one-dimensional integral after the 
introduction of new variables of integration; the final result is 


where 


/.oo 
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I A(Xi - i)\ 3d{ , 


- -L 
4! 


I 3sin3fl 
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> 0 


(5 la) 


(51 A) 


Motivated by the discussion above, we introduce the comparison equation 


dxi 


/.oo 


g - n 


l^g^(AC| - l)d| 


(52a) 


with upstream boundary condition 

g -*■ as X| -► -oo , {S2b) 

for the positive function g - ^X]). This function serves as a lower bound for |/4(X|)| ; 
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I >4(x|)| i 5(^1). The formal solution of (52) is 

00 

Sixi) - • (53n) 

/2-O 

where, for * 0, 1, 2 


a 


n+l 


4!fl 1 

26 (n + 1)(/I + 3/2)6 




/n-0 r«0 


(53A) 


with Oq • 1. It is found numerically that the "tail end" of (53a) is asymptotic to a 
power series whose radius of convergence. can be obtained from the use of the 
ratio test as 

xW - 2.91 - ^ log n . (54) 

Since for x, > 0, (53a) is a series with positive terms, g must become singular at Xj - 
x^i''^ Therefore, | /1(X|)| must become singular at a finite value of x, - x^*^ ^ in view 
of the fact that | /4| ^ g. 

It is easy to show that the solution to the classical Landau equation. 


with upstream condition 


does become singular at 


Q.X I 


A -► e*! as ^1 -o® . 


X| - - 


log J2 . 


(55a) 


(55A) 


(55c) 


Evidently, the averaging of the nonlinear term, as implied by (52a) and (48a). does not 
eliminate the singularity, but merely shifts it downstream. 

Representative numerical results for 0 > ir/4 are presented in figures 7 and 8. It is 
seen from figure 8 that a sectionally linear phase is also consistent with nonlinear 
amplitude equation (44); this observation provides a generalization of trial solution (47). 



- 42 - 


Although the behavior of the magnitude of the amplitude (figure 7) appears 
complicated, an explanation for this behavior is quite simple. It is based on the fact 
that the nonlinear term in (44) depends on the entire history of the amplitude and not 
only on its value at the current location. 

For X, less than about 1.85, equation (48a) for |/1| is applicable since the phase is 
linear. From this equation, we see that the effect of the nonlinear integral is stabilizing 
because sgn(cos 26) <0. As a result, A diminishes to zero (at X| a 1.85) but d/4/dX| is 
nonzero since the latter is balanced by the nonlinear integral that takes into account the 
entire evolution history of the complex amplitude. Hence, A will smoothly pass through 
the origin in the complex plane, thereby suffering a phase jump of ±rr. For x, 
somewhat larger than 1.85, the effect of this phase jump is to render the coefficient of 
d|/l|/dx, in (48a) negative or, equivalently, to make the nonlinear terms destabilizing. 

As X| increases substantially from its value at the first spike (~ 1.85), the effect of 
the initial phase jump begins to show up in the calculation of the nonlinear integral. 
When this occurs to a significant extent, the effect of the phase jump cancels out 
between the cubic terms and the linear terms, and the nonlinear integral is again 
stabilizing. Thereby the second spike (or phase jump) near x, s 3.05 is generated. 
This entire "quasi-periodic" scenario, in which the nonlinear terms are stabilizing or 
destabilizing over certain intervals of the x,-axis. results in the successive spikes of 
figure 7 and in the phase jumps shown in figure 8. The overall trend is a substantial 
increase in the magnitude of A (say. for x, < 4.5), although there are narrow regions in 
the streamwise direction in which the amplitude may be small. This streamwise 
modulation of the disturbance is possible since the nonlinear interaction occurs in the 
form of an integral that is sensitive to upstream values of the amplitude. 

We emphasize that the solutions for the complex amplitude were obtained by 
solving (44) numerically; trial solution (47) and equation (48a) for \A\ provide a 
convenient framework for the interpretation of the numerical results. 
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Finally, when 9 ■ ff/4, the lowest-order nonlinear terms, as given by (42a), vanish. 
In this case, the nonlinear interaction between two oblique modes occurs at higher order 
in the small amplitude parameter. In order to understand this result, it is necessary to 
recall that (in linear theory for a slightly unstable mode) the velocity component of an 
oblique mode normal to the wave front is finite, whereas the component along the wave 
front is singular (in the critical layer). 

When 9 - rr/4, the included angle between the oblique modes is ff/2. so that the 
component of the fluid velocity normal to one family of wave fronts will be completely 
dominated by the tangential velocity along the second family of wave fronts. 
Therefore, the nonlinear interaction between these modes is expected to be weaker than 
in the general case when both waves contribute (roughly equally) to the tangential 
velocity component along the wave fronts. 

To summarize, the nonlinear interaction between two oblique modes in a supersonic 
mixing layer is governed by a nonlinear integro-differential equation of the Hickernell- 
Goldstein type. For each propagation angle, 9 , the magnitude of the scaled amplitude, 
|/4|, is given by a universal curve (in other words, the only essential similarity 
parameter is 9 ). When 9 < rr/4, the disturbance undergoes a monotonic and sudden 
increase in amplitude, which terminates in a singularity at a finite streamwise location. 
On the other hand, for 9 > rr/4, the streamwise development of |/l| is characterized by 
a series of modulations arising from the dual nature of the nonlinearity, which can be 
stabilizing or destabilizing, depending on the streamwise position. 

The authors are grateful for the financial support provided by the NASA Lewis 
Research Center under grant NAG3-78I. One of us (TFB) has benefited from 
stimulating discussions with Drs. M. E. Goldstein, L. Hultgren, and S. S. Lee on several 
aspects of this problem. TTie help of Sang Soo Lee, who provided some of his 
handwritten notes and a computer code, is deeply appreciated. Both items were used to 
provide an independent check on the validity of the results presented herein. 
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Appendix A. Behavior of the solution in the main layers at the edge of the 
critical layer 

The solution in the two main layers (denoted by ±) is given in §4. In order to 
derive the behavior of this solution when {y - y^) is small, we first derive an 
approximation for the integral involving (r*)-'. (see (17i>)). By Taylor series, for x - y - 
y^ small, from (146), 
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and 
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If B vanishes, the generalized inflection point criterion of Lees & Lin (1946) is satisfied. 
We assume that 1/f* is analytic at all y ^ y^. 

The above remarks imply that 
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(A 4) 


is an analytic function for all y and is integrable at j; - ±cx>. The denominator under B 
has been chosen with this objective in mind. Thus, 
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In view of these preparatory remarks, the two-term critical layer expansion of the 
two-term main solution follows from (12), (17). and (A 5): 
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where cl*^ is defined by (lid)- Similarly, 
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where S3* and 0+ “ const, are defined by (23a, b), respectively. Of course, in view of 
these equations, the terms in the braces are independent of the label (±). 

Note that it is possible to give the corresponding results for the cross-stream 
velocity component, and for the perturbation density and temperature. Since these 
results are not needed for the derivation of the amplitude equation, we shall omit them 
for brevity. Observe that u and w in the main layers are more singular as tj -» 0 (i.e. 
y Vc) for oblique modes (I f 0) than for two-dimensional modes (f - 0). This is 
really the main reason why the amplitude scalings (i.e. e - c(ff)) are different for the 
corresponding nonlinear problems. 
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Appendix B. Quadratic nonlinearities 

The equations that govern the evolution of the quadratic nonlinearities (i.e. first 
and zeroth harmonics) in the critical layer are 


^0)|j(l/2) ^1/2) _ _ yO)jj(O) ^ 

(B la) 

^0)^(l/2) _ _ ^O)^(O) ^ 

(B 16) 

g.0/2) gMI/2) ..(1/2) 

di dr) d{ ■ ’ 

(B Ic) 


where and are defined by (196) and (356). respectively. These equations may 
be obtained from (1), (34), and (35). Because the critical layer is a thin layer. 

» 0; therefore, without loss of generality, we have set p(*/2) . 0. Thus, the outer and 
main expansions contain no terms of the form [ff'/ 2 (.)(i/ 2 ) j 

Next, we multiply out the nonlinear forcing terms on the right-hand sides of 
(B la, 6) and write them in Fourier-series representation: 
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2 2 n.tn 
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^—2 
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exp[/(/7{ + mi;)] . 


(B 2a, 6) 


where the and 5T'j are the complex Fourier coefficients. Direct calculation shows 
that - 0; these, in view of (B 16) and null initial conditions at -♦ -oo, 

immediately imply (40a). 

We recall that the lowest order critical-layer solutions (e.g. (39)) are linear in the 
vortex sheet parameter A and, setting A - 0, we recover the corresponding critical-layer 
solutions given by Goldstein & Lee (1992). Thus, for the sake of brevity, we write 
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(GD- 
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^^0 + 9 


dW, 
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dr, 


+ c.c. 


(B 3a) 




where (GL) denotes the terms that are available, at least in principle (see (5.11), (5.13), 
(5.25), and (5.26) in Goldstein & Lee (1992)). Of course, the notation is very different, 
so it is not easy to make a symbol-by-symbol comparison; nevertheless, even a quick 
glance reveals that the solutions are equivalent because they contain the same basic 
information. Note also that c.c. refers to the complex conjugate of the terms that have 
A as a factor in them. 

After using (B 3) and (B 2) in the governing equations, (B 1), and solving these for 
the Fourier coefficients of the velocity components, we find 
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(B 4c) 
(B 4d) 
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Recall that -^i) " 2...) are defined by QAb) and behave as 

for jj -♦ too. Thus, it is possible to match the entire critical-layer solution to that in the 
main layers as given by the linearized theory of §4. Of course, this asymptotic 
matching is carried out only to the order of accuracy of the analysis. Also note that 

when the expressions in (416-A) are multiplied out using (39c) and (B 4), (/fc - 1. 2 

7) become polynomials of degree no greater than two in the vortex sheet displacement 
parameter. A. The evolution equation for the amplitude is expected to depend only on 
even powers of A. since the mirror image of the wavy wall, (6g), cannot change the 
stability characteristics of the oblique modes. Details may be found in Gartside (1995). 

Our final remark concerns the validity of (406). This result follows immediately 
from (B Ic) under the observation that - 0. This comes from the initial 

conditions at - -oo and asymptotic matching. (Note that can be shown to satisfy 
(3/9x| + 2/ai7)(92Ty/^)/9jj2) . 0; we omit the details, however.) 


Appendix C. An illustrative calculation for the phase jump arising from the 
nonlinear terms 

Let us consider the contribution of the nonlinear forcing term, Ty, to the phase- 
jump equation. After using (B 4a) in (416) and solving (37c) for we obtain 

(7) Us?-fi-oL^TQ r ~ 

Q\.\iv. ^i) - 22t7 — ^ ^i) exp[/otiK 2 x, - x,)]dx, 


22i7 — ^ exp(-/a7?x,) 



-^(^i)|^l2(’?. ^i)dx, 


+ A (...) . (C 1) 

where, for the purposes of illustration, we consider only the first two terms in (C 1). 
These are proportional to A^. Recall that /„(i?. x,) (n - 0, 1, ...) are defined by (24c) 
and is the contribution of T7 to 0| | . 
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In order to evaluate 


r 

Gu(>?. Jc,)dij . (C 2) 

J-OO 


we substitute the definition of /2 into (C 1) and use the well-known result 

r OO 


2n 8{x) 
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(C 3) 


where 6(x) is the Dirac delta function. The triple integral in (C 2) is thus reduced to a 
double integral; the latter is further reduced to a single integral by the substitution 
property of the delta function. The final result is 

•^1 

(xi - - X|)dX| 

-oo 

+ A (...) . (C 4) 

to which the second term on the right-hand side of (Cl) has actually made a null 
contribution. We emphasize that this single integral represents a new term that has not 
appeared previously in the literature for the oblique-mode self-interaction problem. 
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